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Abstract 

Random matrix theory (RMT) is based on two assumptions: (1) matrix-element independence, 
and (2) base invariance. Most of the proposed generalizations keep the first assumption and vio- 
late the second. Recently, several authors presented other versions of the theory that keep base 
invariance on the expense of allowing correlations between matrix elements. This is achieved by 
starting from non-extensive entropies rather than the standard Shannon entropy, or following the 
basic prescription of the recently suggested concept of superstatistics. We review these generaliza- 
tions of RMT and illustrate their value by calculating the nearest-neighbor-spacing distributions 
and comparing the results of calculation with experiments and numerical-experiments on systems 
in transition from order to chaos 
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I. INTRODUCTION 



In classical mechanics, integrable Hamiltonian dynamics is characterized by the existence 
of as many conserved quantities as degrees of freedom. Each trajectory in the corresponding 
phase space evolves on an invariant hyper-torus In contrast, chaotic systems are 

ergodic; almost all orbits fill the energy shell in a uniform way. Physical systems with 
integrable and fully chaotic dynamics are exceptional. A typical Hamiltonian system shows 
a mixed phase space in which regions of regular motion and chaotic dynamics coexist . 
Systems of this kind are known as mixed systems. Their dynamical behavior is by no means 
universal. If we perturb an integrable system, most of the periodic orbits on tori with rational 
frequencies disappear. However, some of these orbits persist. Elliptic periodic orbits appear 
surrounded by islands. They correspond to librational motions around these periodic orbits 
and reflect their stability. The Kolmogorov-Arnold (KAM) theorem establishes the stability 
with respect to small perturbations of invariant tori with a sufficiently incommensurate 
frequency vector. When the perturbation increases, numerical simulations show that more 
and more tori are destroyed. For large enough perturbations, there are locally no tori in 
the considered region of phase-space. The break-up of invariant tori leads to a loss of 
stability of the system, to chaos. Different scenaria of transition to chaos in dynamical 
systems have been considered. There are three main scenaria of transition to global chaos 
in finite-dimensional (non-extended) dynamical systems: via the cascade of period-doubling 
bifurcations, the Lorenz system-like transition via Hopf and Shil'nikov bifurcations, and the 

nnfi 

transition to chaos via intermittences [3|, U, |5|. It is natural to expect that there could be 
other (presumably many more) such scenaria in extended (infinite-dimensional) dynamical 
systems. 

In quantum mechanics, the specification of a wave function is always related to a certain 
basis. In integrable systems eigenbasis of the Hamiltonian is known in principle. In this 
basis, each eigenfunction has just one component that obviously indicates the absence of 
complexity. In the nearly ordered regime, mixing of quantum states belonging to adjacent 
levels can be ignored and the energy levels are uncorrelated. The level-spacing distribution 
function obeys the Poissonian, exp(— s), where s is the energy spacing between adjacent 
levels normalized by the mean level spacing D. On the other hand, the eigenfunctions a 
Hamiltonian with a chaotic classical limit is unknown in principle. In other words, there is 
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no special basis to express the eigenstates of a chaotic system. If we try to express the wave 
functions of a chaotic system in terms of a given basis, their components become on average 
uniformly distributed over the whole basis. They are also extended in all other bases. For 
example, Berry [6] conjectured that the wavefunctions of chaotic quantum systems can be 
represented as a formal sum over elementary solutions of the Laplace equation in which real 
and imaginary parts of coefficients are independent identically-distributed Gaussian random 
variables with zero mean and variance computed from the normalization. Bohigas et al. 

put forward a conjecture (strongly supported by accumulated numerical evidence) that 
the spectral statistics of chaotic systems follow random-matrix theory (RMT) j^, . This 
theory models a chaotic system by an ensemble of random Hamiltonian matrices H that be- 
long to one of the three universal classes, orthogonal, unitary and symplectic. The theory is 
based on two main assumptions: the matrix elements are independent identically-distributed 
random variables, and their distribution is invariant under unitary transformations. These 
lead to a Gaussian probability density distribution for the matrix elements. The Gaussian 
distribution is also obtained by maximizing the Shannon entropy under constraints of nor- 
malization and existence of the expectation value of Tr (H^H) , where Tr denotes the trace 



and stands for the Hermitian conjugate of H [8|, [lOj. The statistical information about 
the eigenvalues and/or eigenvectors of the matrix can be obtained by integrating out all the 
undesired variables from distribution of the matrix elements. This theory predicts a univer- 
sal form of the spectral correlation functions determined solely by some global symmetries 
of the system (time- reversal invariance and value of the spin) . Time-reversal- invariant quan- 
tum system are represented by a Gaussian orthogonal ensemble (GOE) of random matrices 
when the system has rotational symmetry and by a Gaussian symplectic ensemble (GSE) 
otherwise. Chaotic systems without time reversal invariance are represented by the Gaus- 
sian unitary ensemble (GUE). The dimension (3 of the underlying parameter space is used 
to label these three ensembles: for GOE, GUE and GSE, (3 takes the values 1, 2 and 4, re- 
spectively. Among several measures representing spectral correlations, the nearest-neighbor 
level-spacing distribution function p{s) has been extensively studied so far. According to 
the random matrix theory, the level spacing distribution function in the chaotic phase is 
approximated by the Wigner-Dyson distribution, namely. 
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The coefficients a^? and bjs are determined by the normahzation conditions Pp{s)ds = 
sP,3{s)ds = 1, as ai = n/2, as = 32/7r2, = 2i7367r3, bi = n/A, 62 = 4/7r,and 64 = 
64/97r. For s <^ 1, the distribution function is proportional to s^, which imphes that 
adjacent energy levels cannot approach each other indefinitely because of mixing between 
two extended states. 

The assumptions that lead to RMT do not apply for mixed systems. The Hamiltonian 
of a typical mixed system can be described as a random matrix with some (or all) of its 
elements as randomly distributed. Here the distributions of various matrix elements need 
not be same, may or may not be correlated and some of them can be non-random too. 
This is a difficult route to follow. So far in the literature, there is no rigorous statistical 
description for the transition from integrability to chaos. There have been several proposals 
for phenomenological random matrix theories that interpolate between the Wigner-Dyson 
RMT and banded RM with the (almost) Poissonian level statistics. The standard route of 
the derivation is to sacrifice basis invariance but keep matrix-element independence. The first 
work in this direction is due to Rosenzweig and Porter Hj] . They model the Hamiltonian of 
the mixed system by a superposition of two matrices: a diagonal matrix of random elements 
having the same variance and a matrix drawn from a GOE. Therefore, the variances of the 
diagonal elements total Hamiltonian are different from those of the off-diagonal ones, unlike 
the standard GOE Hamiltonian in which the variances of diagonal elements are twice of the 
off-diagonal ones. Hussein and Sato [l2| used the maximum entropy principle to construct 
such ensembles by imposing additional constraints. Ensembles of band random matrices 
whose entries are equal to zero outside a band of width b along the principal diagonal have 
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16|,ll7|. 



also been often used to model mixed systems 

Another route for generalizing RMT is to conserve base invariance but allow for corre- 
lation of matrix elements. This has been achieved by maximizing non-extensive entropies 



subject to the constraint of fixed expectation value of Tr(^H^H) 18 
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24|. 



Recently, an equivalent approach is presented in 25|, l26|], which is based on the method of 
superstatistics (statistics of a statistics) proposed by Beck and Cohen 27|]. This formalism 
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Deen elaborated and a 
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ied successfully to a wide variety of physical problems, e.g., in 
35|. In thermostatics, superstatistics arises as weighted averages 
of ordinary statistics (the Boltzmann factor) due to fiuctuations of one or more intensive 
parameter (e.g. the inverse temperature). Its application to RMT assumes the spectrum of a 
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mixed system is made up of many smaller cells that are temporarily in a chaotic phase. Each 
cell is large enough to obey the statistical requirements of RMT but has a different distribu- 
tion parameter t] associated with it, according to a probability density f{f]). Consequently, 
the superstatistical random-matrix ensemble describes the mixed system as a mixture of 
Gaussian ensembles with a statistical weight /(r^). Therefore one can evaluate any statistic 
for the superstatistical ensemble by simply integrating the corresponding statistic for the 
conventional Gaussian ensemble. 



II. NONEXTENSIVE GENERALIZATION OF RMT 



In 1957 Jaynes [39|] proposed a rule, based on information theory, to provide a constructive 
criterion for setting up probability distributions on the basis of partial knowledge. This leads 
to a type of statistical inference which is called the maximum-entropy principle (MaxEnt). 
It is the least biased estimate possible on the given information. Jaynes showed in particu- 
lar how his rule, when applied to statistical mechanics, leads to the usual Gibbs' canonical 
distribution. The core of the MaxEnt method resides in interpreting entropy, through the 
Shannon axioms, as a measure of the "amount of uncertainty" or of the "amount of infor- 
mation that is missing" in a probability distribution. This was an important step forward 
because it extended the applicability of the notion of entropy far beyond its original roots in 
thermodynamics. Balian considered the application of MaxEnt to the random-matrix the- 
ory by maximizing the Shannon entropy under constraints of normalization and existence 
of the expectation value of Tr(H"''H). In this section, we consider possible generaliztion of 
RMT by extremizing two different entropies, namely Tsallis' and Kaniadakis', rather than 
Shannon's entropy. The extremization is again subject to the constraint of normalization 
and existence of the expectation value of Tr(H^H). For completeness, we start by a brief 
review of the conventional random-matrix theory. 



A. RMT from Shannon's entropy 



Balian [10|] derived the weight functions -P(H) for the random-matrix ensembles from 
MaxEnt postulating the existence of a second moment of the Hamiltonian. He applied the 
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conventional Shannon definition for tlie entropy to ensembles of random matrices as 

Ssci. = -j dHP(H)lnP(H) (2) 

and maximized it under the constraints of normalization of -P(H) and fixed mean value of 
Tr(H"''H). The latter constraint ensures basis independence, which is a property of the trace 
of a matrix. Then, the distribution P{H) is determined from the extremum of the functional 

i^sch = ^sch -^JdH P(H) -r]j rfHP(H)Tr {H^H) , (3) 

where ^ and t] are Lagrange multipliers. Its maximum is obtained equating its functional 
derivative to zero. He obtained 

(H) = ^ exp [-r/Tr {U^U) ] , (4) 

where Z{ri) = J exp [— r/Tr (H^H)] dH. 

It is easy to see that the joint distribution of matrix elements obtained in Eq. (4) satisfies 
the two conditions of RMT, namely uncorrelated matrix-elements and Base independence. 
The first condition follow since the distribution (4) is a Gaussian distribution with inverse 
variance 1/21], since Tr(H''"H) = J2iLi (j^ii''^ + '^J2y=oJ2i>j {^ij^^ ■> where all the four 
matrices H^'^^ with 7 = 0,1,2,3 are real. This allows the factorization of P/? (H) into 
products of terms depending only on the individual matrix elements. Therefore, the matrix 
elements of H are independent. Base independence follows from the fact that the distribution 
(18) depends on H in the combination Tr(H^H). Indeed, if two matrices A and B that 
express the same operator in two different bases are related by a similarity transformation 
B = T^^AT, then such operators have the same trace. 

The joint distribution of eigenvalues Ei immediately follows from Eq. (3). With H = 
U^^EU, where U is the global unitary group and E = diag(i?i, ■ ■ ■ , En) the volume element 
da has the form 

<m = \ An {E)f dEdfiiV), (5) 

where AAr(E) = Y[n>m(^n — Em) is the Vandermonde determinant and dfi(\J) the in- 
variant Haar measure of the unitary group j^, ^. Integrating over U and noting that 
Tr (H"!"!!) =TrE^ yields the joint probability density of eigenvalues in the form 

Pf,{v, E,,---, En) = Cf,l[ \En - Emf exp ( -r^ Ef] , (6) 



n>m \ i=l 



where is a normalization constant. All of the spectral properties of the Gaussian random- 
matrix ensemble can be obtained from Eq. (6). However, this is not an easy task. Lacking 
simple exact results, and guided by the case N = 2, Wigner proposed a form for the nearest 
neighbor spacing (NNS) distribution p{s) of eigenvalues. This "Wigner surmise" , originally 
stated for (3 = 1, has the form 



Pfiis,v) = T^uf^, \ /oi (Tr)^exp f-^ ) . (7) 



r[(/3+l)/2] V2y 2 

The parameter r) is determined by the condition of unit mean spacing, s pfi{s)ds = 1, as 

^ ^ r^l(/3 + 2)/2| 
' P|(/3+l)/2|' ^' 

Although the Wigner surmise is strictly valid for two-dimensional ensembles, it is an accurate 
approximation for ensemble with arbitrarily large A^. To demonstrate the accuracy, we 
expand this distribution for the case of /3 = 2 in powers of s to obtain 



P2(s) = ( 1 - + . . . ) ^ 3.242s' - 4.128s' + ■ ■ ■ , (9) 



while thepower-series expansion of the corresponding exact distribution for ensembles with 
A^ oo y yields 

P2,exact(s) = " + ■■■ = 3.2905' - 4.329s' + ■ • • . (lO) 

The Wigner surmise has been successfully applied to the NNS distributions for numerous 
chaotic systems. 

B. RMT from Tsallis' entropy 

The past decade has witnessed a considerable interest devoted to non-conventional sta- 
tistical mechanics. Much work in this direction followed the line initiated by Tsallis' seminal 



paper 40|. The standard statistical mechanics is based on the Shannon entropy measure 
S = —HiPilnpi (we use Boltzmann's constant ks = 1), where {pi} denotes the probabili- 
ties of the microscopic configurations. This entropy is extensive. For a composite system 
A + B, constituted of two independent subsystems A and B such that the probability 
p{A + B)= p{A)p{B), the entropy of the total S{A + B) = S{A) + S{B). Tsallis proposed 
a non-extensive generalization: Sg = {1 — T^ipf) /(g — 1). The entropic index q characterizes 



the degree of extensivity of the system. The entropy of the composite system A + B, the 
Tsalhs' measure verifies 



S,{A + B) = S,{A) + S,{B) + (1 - q)S,{A)S,{B), (11) 

from which the denunciation non-extensive comes. Therefore, Sq{A + B) < Sq{A) + Sg{B) 
if g > 1. This case is called sub-extensive. If g < 1, the system is in the super-extensive 
regime. The standard statistical mechanics recovered for q = 1. Applications of the Tsallis 



formalism covered a wide class of phenomena; for a review please see, e.g. 41 1. 

The Tsallis entropy is defined for the joint matrix-element probability density P^(g, H) 

by 

S, [PtsAq, H)] ={l- j dH [PtsA^I. H)]^) I {q-l). (12) 

We shall refer to the corresponding ensembles as the Tsallis orthogonal ensemble (TsOE), the 
Tsallis Unitary ensemble (TsUE), and the Tsallis symplectic ensemble (TsSE). For q 1, 
Sq tends to Shannon's entropy, which yields the canonical Gaussian orthogonal, unitary or 
symplectic ensembles (GOE, GUE, GSE) |8|,Q. 

There are more than one formulation of non-extensive statistics which mainly differ in 
the definition of the averaging. Some of them are discussed in 42|. We apply the most 



recent formulation [43|]. The probability distribution Pts,/3('?, H) is obtained by maximizing 
the entropy under two conditions, where aj3 is a constant. The optimization of Sq with these 
constraints yields a power-law type for 

Pp{q, H)P^,Aq, H) = [1 + (g - l)vq {Tr (RtR) - a^}] , (13) 

where rjq > is related to the Lagrange multiplier r] associated with the constraint of fixed 
Tr(HtH) by 

V, = V/ J ^^H[PTs,/5(g,H)f , (14) 

and 

Zq = j <m[l + {q- l)?jq {Tr (HtH) - cr^}]"^ . (15) 

It turns out that the distribution (13) can be written hiding the presence of in a more 
convenient form 

PTs,/3(g, H) = Z-^ [l + {q- l)r/,Tr (fitH)] , (16) 
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where 



and 



Vq 



JdU [P^(g,H)r + (l-gW' 



(17) 



rfH [1 + (g - (H^H)] . (18) 

The probabihty density Pts,/3(q', H) depends on H through Tr(H^H) and is therefore invari- 
ant under arbitrary rotations in the matrix space. This ensures base invariance. It decays 
by a power law as the square of any matrix element tends to infinity in contrast with the 
Gaussian decay of the distribution function of the conventional random matrix ensembles. 

We now calculate the joint probability density for the eigenvalues of the Hamiltonian H. 
Expressing the volume element in the matrix-element space in the form (5) and integrating 
over the "angular" variables, on obtains 



N 



G 



n 



Ts,/3 I I 
n>m 



N 



i=l 



(19) 



where Cts,/3 is a normalization constant. 

In order to obtain a generalization of Wigner's surmise, we consider the special case of 
= 2. In this case. 



1 + (g - l)r/, ( 2e' + -s' 



9-1 



(20) 



where e = {Ei + E2)/2, s = \Ei — E2\. For this case, the distribution fl20|) has to be 
complemented by the auxiliary condition that the quantity inside the square bracket has to 
be positive. We here consider the case of g > 1 where no limitations are imposed on the 
values of the variables e and s, and refer the reader interested in the other case of g < 1 to 
Ref. [21]. The NNS distribution is obtained by integrating fl^U]) over e from —00 to 00. 



PTsA^, s) = aTs,/3S^ [1 + 6ts,/3S^] ^ 



q-l ' 2 



(21) 



where aTs,/3 is a normalization coefficient and &ts,/3 is obtained for the requirement of unit 
mean spacing. Explicitly, 



0'Ts,l3 



^,(/3+l)/2 p /j^ 1 



r m r 



and 6ts,/: 



2 



(22) 
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The second moment of the distribution (s^) = s^pp{s)ds is given by 




It diverges unless q < q^o = 1 + 2/ (/3 + 4), which is equal to 1.40, 1.33 and 1.25 for the 
orthogonal, unitary and symplectic ensemble, respectively. This imposes physical bound 
on the admissible values of g, because (s^) has to be finite in order to force condition that 
Tr(H"''H) has to be finite. At higher values of the entropic index, non-extensive statistics 
does not apply to the random matrix model. The peak of the distribution in Eq. (1211) is 
located at sp = \Jl3/h'Xs,i3 [—1 — /9 + 1/(1 — Q')] • It moves from si = 0.798 to si = 0.368, 
from $2 = 0.886 to $2 = 0.408 and from S4 = 0.940 to S4 = 0.671 as g = 1 (the standard 
Wigner's surmise) to goo- Neither reaches 0, the peak position of the Poisson distribution 
exp(— s) of the integrable systems. The proposed non-extensive ensemble in the three cases 
of /3 = 1, 2 and 4 evolve the shape predicted by the corresponding Wigner surmise towards 
the Poisson distribution, but never reach it. 



C. RMT from Kaniadakis' entropy 

In this section, we consider a possible generalization of RMT based on an extremization 



of Kaniadakis' K-entropy |44j . l45l . l46l ]. This entropy shares the same symmetry group of 
the relativistic momentum transformation and has applications in cosmic-ray and plasma 
physics. For the matrix-element probability distribution function, it reads 

[k, P^^p{K, Yi)] = -^ f dH [Pk,/3(«:, m'^^ - [^K,/3(«:, H)]^"'^) (24) 

Ik J K V — K J 

with K a parameter with value between and 1; the case of k = corresponds to the 
Schannon entropy. Here, a is a real positive parameter. Kaniadakis has considered two 
choices of a, namely a = 1 and a = Z, where Z is the generalized partition function. We 
here adopt the second choice. The matrix-element distribution Pk,/3('^;H) is obtained by 
extremizing the functional 

FK = S^-VKj dH PkA^, H)Tr {H^H.) , (25) 

where t^k is a Lagrange multiplier. One arrives to the following distribution 

PkA^^, H) = ^ exp^,} [-r^KTr (Htfi)] , (26) 
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where 



Z^ = Jdn exp^^y [-r/xTr (HtH)] . 
Here exp|^| [x] is the K-deformed exponential ^] which is defined by 

expj^l [x] = ^Vl + K^x^ + Kxj = exp ^— arcsinh kx^ . 
Concerning the asymptotic behavior of PK,f3{K, H) we easily obtain that 

Pk,/3(«:,H)~ |«r?KTr(HtH)| 



(27) 



(28) 



(29) 



as the square of any of the matrix elements tends to infinity 

The probability density Pk(k, H) depends on H through Tr(H^^H) and is therefore in- 
variant under arbitrary rotations in the matrix space. Using Eq. (26) and integrating over 
the "angular variables" U yields the joint probability density of eigenvalues in the form 



Pk,/3{k;Ei,...,En) = Ck,i3 Y\_ \En - Emf exp^^ 



n>m 



N 
i=l 



(30) 



where Ck,/3 is a normalization constant. 

In order to obtain a generalization of the Wigner surmise, we consider the case of two- 
dimensional random-matrix ensemble where N = 2 and n = 3 and restrict our consideration 
to the orthogonal ensemble with /? = 1. In this case, Eq. (!30l) reads 



Pk,/3(/«; £, s) = ^ ^^^t \ (f^Vv^^^^ 



exp{«} 



-r/K ( 2e^ + -s^ 



(31) 



^ - I'D 

where e = [Ei + E2) /2, s = \Ei — E2\, and B{a,b) = r(a)r(6)/r(a + 6) is the Beta function 



471]. The NNS distribution is obtained by integrating (l3Ti) over e from —00 to 00. This can 



be done by changing the variable e into x = exp [—■iarcsinh(K?7K ■5^/2)], integrating by parts, 
and then replacing the variable x by another new variable, y = exp(/ta;). The resulting 
integral can be solved by using the following identity 



y-\y + (y - dy = u^+'^-'b (a - /i - z., /.) 



-u, X — jj — u ; \ — u; — — I , (32) 

uj 
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for \l3/u\ < 1 and < ;U < A — z/, where 2Fi{i',fi; \;x) is the hypergeometric function. Thus, 
after straightforward calculations we can express the NNS distribution as 

PkA'^,s) = -2(1 + ^^,,W2-iM^rcsinH^nKsV2) B - i I 

(--, ---;- + !; _e-2arc.mh(.„KsV2)^ , (33) 

y 2 /t^ 2 A€ J 

The condition of unit mean spacing defines the quantity r^x as 



Vk 



(34) 



We note that the function _B (^ — |, |) diverges at k = = 1/2, which serves as an upper 
limit for the range of variation of k. We also note that the mean square spacing diverges at 
Kj —— i^oo ~~ 2y'5. 

III. SUPERSTATISTICAL GENERALIZATION OF RMT 

Let us first recall the basic idea underlying superstatistics. We will then proceed to 
construct a generalization of RMT in the spirit of superstatistics. 

A. Beck and Cohen's superstatistics 

Consider a complex system in a nonequilibrium stationary state. Such a system will 
be, in general, inhomogeneous in both space and time. Effectively, it may be thought to 
consist of many spatial cells, in each of which there may be a different value of some relevant 
intensive parameter, e.g. the inverse temperature /3. Beck and cohen assumed that this 
quantity fluctuates adiabatically slowly, namely that the time scale is much larger than the 
relaxation time for reaching local equilibrium. In that case, the distribution function of the 
non-equilibrium system consists in Boltzmann factors exp {—(3H) that are averaged over the 
various fluctuating inverse temperatures 

F(H) = jgipf-^^^ip (35) 

where z{l3) is a normalizing constant, and g{[3) is the probability distribution of (5. Let us 
stress that (3~^ is a local variance parameter of a suitable observable, the Hamiltonian of 
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the complex system in this case. Ordinary statistical mechanics are recovered in the limit 
g{l3) 6 {j3). In contrast, different choices for the statistics of may lead to a large variety of 
probability distributions F{H). Several forms for g{f3) have been studied in the literature, 



e.g. 



27 



32 



48| . In particularly, Tsallis statistics is generated from Eq. fl55]) if /3 is a 



chi-squared random variable 



1 / \ 1^/2 

1 V ^ ' 



where /3o is the mean value. A chi-squared distribution derives from the summation of 
squares of v Gaussian random variables Xi, . . . , Xy due to various relevant degrees of 
freedom in the system, where the Xj are independent, and (Xj ) = 0. If (5"^ rather than /? 
is the sum of several squared Gaussian random variables, the resulting distribution gij]) is 
the inverse distribution given by 



where again (3q is the average of (3. Instead of being a sum of many contributions, the random 
variable (3 may be generated by multiplicative random processes. Then In /3 = Yl'i=i 1^ -^i 
a sum of Gaussian random variables. Thus it is log-normally distributed, 

^(/5) = ^i-e-[^°(^/-)lV-^ (38) 



which has an average fi^/w and variance fi'^w{w — 1), where w = exp(f^). 



B. RMT within superstatistics 

To apply the concept of superstatistics to RMT, one assumes the spectrum of a (mixed) 
system as made up of many smaller cells that are temporarily in a chaotic phase. Each 
cell is large enough to obey the statistical requirements of RMT but is associated with 
a different distribution of the parameter rj in Eq. (4) according to a probability density 
/(r/). Consequently, the superstatistical random-matrix ensemble used for the description 
of a mixed system consists of a superposition of Gaussian ensembles. Its joint probability 
density distribution of the matrix elements is obtained by integrating the distribution given 
in Eq. (4) over all positive values of t] with a statistical weight f{ri), 

(39) 
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Despite the fact that it is hard to make this picture rigorous, there is indeed a representation 



which comes close to this idea 



50|. 



Beck, Cohen and Swinney [48|] proposed the derivation of superstatistics starting from 
time-series. The idea is that superstatistical thermostatics resuhs as a convolution of two 
statistics, one characterized by the Boltzmann factor and the other corresponding to inverse- 
temperature fluctuations. This requires the existence of two relaxation times. A justification 
for the use of the above-mentioned superstatistical generalization of RMT in the study of 
mixed systems, is given in 51]. It is based on the representation of their energy spectra 
in the form of discrete time series in which the level order plays the role of time. The 
representation of the suitably transformed eigenvalues of a quantum system as a time series 



las recently a. 



52 



53 



54 



55 



lowed to determine the degree of chaoticity of the dynamics of the system 



561] . Reference 51i considers two billiards with mushroom-shaped boundaries 



as representatives of systems with mixed regular-chaotic dynamics and three with the shape 
of Limagon billiards, one of them of chaotic and two of mixed dynamics. The quantum 
eigenvalues and statistical properties of the eigenfunctions were obtained experimentally 
by exploiting the equivalence of the Schrodinger equation of a plane quantum billiard and 
the Helmholtz equation for the electric field strength in a cylindrical microwave resonator 
for wave lengths longer than twice the height of the resonator. The billiards with mixed 
dynamics have classical phase spaces of different structures for the two families of billiards. 
The "time-series" analysis of their spectra indeed manifests the existence of two relaxation 
lengths. The short one, which is defined as the average length over which energy fluctuations 
are correlated, is of the order of the mean level spacing. The long one, which is by an order 
of magnitude larger, characterizes the typical linear size of the heterogeneous domains of 
the total spectrum. 

The new framework of RMT provided by superstatistics should now be clear. The local 
mean spacing is no longer uniformly set to unity but allowed to take different (random) values 
at different parts of the spectrum. The parameter t] is no longer a fixed parameter but it 
is a stochastic variable with probability distribution f{r]). The observed mean level spacing 
is just the expectation value of the local ones. The fluctuation of the local mean spacing is 
due to the correlation of the matrix elements, which disappears for chaotic systems. In the 
absence of these fluctuations, /(r/) = 6{ri — r/o) and we obtain the standard RMT. Within 
the superstatistics framework, we can express any statistic cxlE) of a mixed system that can 
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in principle be obtained from the joint eigenvalue distribution by integration over some of 
the eigenvalues, in terms of the corresponding statistic a^^\E,r]) for a Gaussian random 
ensemble. The superstatistical generalization is given by 

POO 

a{E)= / /(r7)a(^)(i?,r/)dr7. (40) 
Jo 

The remaining task of superstatistics is the computation of the distribution /(r^), which has 
been introduced in Eq. (|39l) . The time series analysis in Ref. SlJ allows to derive the sape 
the parameter distribution f{ri). The obtained distribution agrees better with the inverse 
distribution given by Eq. (I57|) rather than the other two distributions (1361) and (1381) . We 
have already mentioned that the distribution of the superstatistical parameter rj yields 
Tsallis statistics for RMT, which is considered in the previous section. The log-normal 
distribution does not lead to simple analytical results for the important level statistics like 
the NNS distribution. For these reasons, we shall confine our further consideration to the 
case of inverse distributed superstatistical parameter rj. 



C. Superstatistical generalization of Wigner's surmise 



The superstatistics generalization of NNS distribution for an ensemble belonging to a 
given symmetry class is obtained by substituting the NNS distribution of the corresponding 
Gaussian ensemble pp{s,r]) for a^^\E,r]) in (140|) and integrating over the local mean level 
spacing rj 

/■oo 

PssAs)= f{v)pp{s,v)dv- (41) 
Jo 

For an inverse distribution of r], given by Eq. fl37j) . one obtains the following superstatis- 
tical NNS distribution 



(42) 



r(|)r(i±2) 

where Km{x) is a modified Bessel function [47|] and r]Q again is determined by the requirement 
that the mean-level spacing (s) equals unity, 

>(^)r(i + f)^^ 



IGtt 



^0 



(43) 



r(|)r(iM) 

The inverse distribution of rj follows when the quantity r]~^ is the sum of u squared 
Gaussian random variables. If we take this assumption literarily, we must restrict u to take 
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positive integer values. As the transition from integrability to chaos is known to proceed 
continuously, we have to relax this condition and allow v to take any real value between 1 
and oo. Let us restrict our following consideration to the case of orthogonal symmetry with 
/? = 1. Using the asymptotic expression of the modified Bessel function [47,], we easily find 

limpss,i(i^,s) = ^se--^'/', (44) 

which is the Wigner surmise, as required. The other limit of u ^ 1 yields the semi-Poisson 
distribution 

PSemiPoisson(s) = 4se"^^ (45) 

which is known to provide a satisfactory description for the spectra of pseudointegrable 



systems such as planar polygonal billiards, when all their angles are rational with tt 



57j. The 



motion of the corresponding classical systems in phase space is not restricted to a torus like 



for integrable systems, but to a surface with a more complicated topology [58|. We therefore 
conclude that the assumption that the inverse square of the variance of matrix elements 
as an inverse variable allows one to model the transition out of chaos (corresponding to 
z/ ^ 1) until the system reaches the state of quasi- integrability as the effective number u of 
components of t]~^ approaches 1. If one allows u to take lower values, then the distribution 
l2l) tends to the Poisson distribution as ^ — 1; 

Pss,ii-l,s)=e-\ (46) 



We there conclude that formula (1421) can provide a model for describing the stochastic 
transition all the way from integrability to chaos passing by the stage of quasi-integrability. 



IV. COMPARISON WITH NUMERICAL EXPERIMENT 

The NNS distributions Pts,/3(q', s) and PK,/3ii^, s) obtained above when the entropy is given 
by the Tsallis and Kaniadakis entropies, respectively, as well as the superstatistical distribu- 
tion pss,i3{^, s) describe evolution of the spacing distribution from the Wigner shape to the 



Poissonian. They can be useful for describing systems with mixed regular-chaotic at 



east 



when more familiar distributions such as Berry and Robnik's or Brody's distribution |611.l62| 
do not fit the data satisfactorily. We shall demonstrate this by using these distributions to 
fit the NNS distribution of levels of a number of mixed systems. 
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A. Mushroom billiards 



Billiards can be used as simple models in the study of Hamiltonian systems. They consist 
of a point particle which is confined to a container of some shape and reflected elastically 
on impact with the boundary. The shape determines whether the dynamics inside the 
billiard is regular, chaotic or mixed. The best-known examples of chaotic billiards are the 
Sinai billiard (a square table with a circular barrier at its center) and the Bunimovich 

n 

stadium (a rectangle with two circular caps) [59|. Neighboring parallel orbits diverge when 
they collide with dispersing components of the billiard boundary. In mixed billiards, some 
neighboring parallel orbits converge at first, but divergence prevails over convergence on 
average. Divergence and convergence are balanced in integrable billiards such as circles and 
ellipses. 

Recently Bunimovich introduced the so-called 'mushroom' billiard [gO] with the novel 
feature of a well-understood divided phase-space comprising a single integrable region and 
a single ergodic one. We restrict ourselves here to mushroom billiards which consist of a 
semicircular region, the 'hat' and a 'stem', which is symmetrically attached to its base. As 
the width of the stem varies from zero to the diameter of the hat, there is a continuous 
transition from integrability (the semicircle billiard) to ergodicity (in case of a rectangular 
stem the stadium billiard). In mushroom billiards, the regular region has a well-defined 
semicircular border. It is composed of those trajectories in the hat that never cross this 
border and therefore remain in the hat forever. Their integrability is due to the conservation 
of the reflection angle for collisions with the semicircular boundary. The chaotic component 
consists of trajectories that enter the stem of the mushroom billiard. Two mushroom billiards 
have been recently investigated experimentally 63|]. The ratio of the width of the stem to 



the diameter of the hat is 1:3 for the small mushroom billiard and 2:3 for the large. Both 
billiards have mixed dynamics with classical phase spaces of different structures for the 
two billiards. The degree of chaos, which is the measure of all chaotic parts of the phase 
space, is 45.5 % and 82.9 %, respectively. Both systems have been studied experimentally 
in the quantum limit exploiting the analogy between a quantum billiard and a fiat cylindric 
microwave billiard of the same shape. The electromagnetic resonances in each microwave 
cavity can directly be associated with quantum states in the corresponding quantum billiard. 
For the evaluation of statistical measures, a sufficiently large number of resonances is needed. 
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The first 780 resonances could be detected in the small billiards and 938 in the large one. 
The quantum eigenvalues were obtained experimentally by exploiting the equivalence of 
the Schrodinger equation of a plane quantum billiard and the Helmholtz equation for the 
electric field strength in a cylindrical microwave resonator for wave lengths longer than twice 
the height of the resonator. To compare the statistical properties of the eigenvalues with 
universal predictions considered in the present paper, they have to be rescaled to unit mean 



spacing. This is done by an unfolding procedure using Weyl's formula [6J] , which relates the 
billiard area and circumference to the number of resonance frequencies below a given one. 

We compared the resulting NNS distributions given in Eqs. ( 12T1) . ( l33i) and ( H2l) with the 
experimental ones for the mushroom billiards and the two Limagon billiards with mixed 
dynamics. In Fig. 1 the experimental results for pi{s) are shown by histogram together 
while the distributions obtained by starting with the Tsallis and Kaniadakis entropies and 
the superstatistical distributions are shown by the dashed, dashed-dotted and solid lines, 
respectively. The best fit values of the parameters are given in Table I, together with the 

values calculated by 

i=l 

where is the number of experimental spacings and X stands for Ts (Tsallis), K (Kani- 
adakis) or SS (superstatistics). 

B. Random binary networks 

As another example of mixed systems, we consider a numerical experiment by Gu et al. 
651] on a random binary network. Impurity bonds are employed to replace the bonds in an 
otherwise homogeneous network. In such a network, there exist a lot of geometric resonances 
randomly distributing in the resonant area. Based on the Green's-function formalism, the 
eigenvalues of Green's-matrix are solved, the sequence of which forms the resonance spectrum 



The authors of Ref. 65| numerically calculated more than 700 resonances for each sample. 
For each impurity concentration p, they considered 1000 samples with totally more than 
700 000 levels computed. Their results for four values of concentration p are compared with 
both the Tsallis, Kaniadakis and superstatistical NNS distributions in Fig. 2. The best fit 
values of the parameters as well as the corresponding value are given in Table II. The high 



18 



statistical significance of the data allows us to assume the advantage of the superstatistical 
distributions for describing the results of this experiment, as compared to the other two 
distribution families. 

V. CONCLUSION 

Random matrix theory is the statistical theory of random matrices, whose entries fluctu- 
ate as independent Gaussian random numbers. The matrix-element distribution is obtained 
by extremizing Shannon's entropy subject to the constraint of normalization and constant 
trace of the square of the matrix. The latter constraint renders the matrix-element dis- 
tribution base independent. While most of the previously proposed generaliztion of RMT 
violate base invariance, the ones reviewed in this paper conserve it. Non-extensive gen- 
eralizations extremize nonextensive entropies such as Tsallis' or Kaniadakis', rather than 
Shannon's. Superstatistical generalizations, on the other hand, allow the fluctuation of the 
mean local density of states which is flxed in the standard theory. These generalizations of 
RMT, seen from the present perspective, may bear interest per se because of the additional 
nontrivial fluctuations introduced in a simple model. In addition, they may constitute a use- 
ful statistical paradigm for the analysis of the spectral fluctuations of systems with mixed 
regular-chaotic dynamics. For this purpose, simple analytical expressions are derived in each 
case for the nearest neighbor level distributions, being among the most popular characteris- 
tics of spectral fluctuation. The formalism has been checked by the analysis of experimental 
resonance spectra of mixed microwave billiards and geometrical resonances in random bi- 
nary networks. The predictions of the three models satisfactorily describe the experimental 
trends of the evolution of NNS distributions during the transition out of chaos. The con- 
sidered experimental data agree better in most cases with the corresponding distributions 
predicted by the superstatistical approach when the fluctuating parameter has an inverse 

distribution. 
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TABLE I. Best-fit parameters for the experimental NNS distribution of resonances in the 

small and large mushroom billiards. The corresponding ^ values are given in parentheses. 
Distribution Small billiard Large billiard 

Tsallis q = 1.336 (0.0189) q = 1.221 (0.0031) 

Kaniadakis k = 0.423 (0.0159) k = 0.017 (0.0069) 
Superstatistical i/ = -0.441 (0.0026) = 2.31 (0.0018) 

TABLE IL Best-fit parameters for the numericcal-experimental NNS distribution of ge- 
ometrical resonances in the binary random network with different impurity concentrations 

p. The corresponding values are given in parentheses. 

Distribution p = 0.1 p = 0.2 p = 0.3 p = 0.4 

Tsallis q = 1.380 (0.0034) q = 1.322 (0.0022) q = 1.263 (0.0015) q = 1.219 (0.0009) 

Kaniadakis k = 0.444 (0.0232) k = 0.421 (0.0067) k = 0.398 (0.0059) k = 0.012 (0.0051) 

Superstatistical u = -0.188 (0.0021) u = 0.617 (0.0002) u = 1.76 (0.0002) u = 3.12 (0.0003) 
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FIG.l NNS distributions of resonances in mushroom billiards (histograms), measure by 
Friedrich et al. [gsI compared with the Tsallis (solid lines), Kaniadakis (dashed) and super- 
statistical (dashed-dotted) NNS distributions. The best-fit parameters are given in Table 
1. 
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FIG. 2 NNS distributions of geometrical resonances in random network (dots), calculated 
by Gu et al. |65| compared with the Tsallis (solid lines), Kaniadakis (dashed) and super- 
statistical (dashed-dotted) NNS distributions. The best-fit parameters are given in Table 
11. 
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